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A method to determine the quantum state of a scalar held after 0(4)-symmetric bubble nucleation 
has been developed recently. The method has an advantage that it concisely gives us a clear picture of 
the resultant quantum state. In particular, one may interpret the excitations as a particle creation 
phenomenon just as in the case of particle creation in curved spacetime. As an application, we 
investigate in detail the spectrum of quantum excitations of the tunneling held when it undergoes 
false vacuum decay. We consider a tunneling potential which is piece-wise quadratic, hence is simple 
enough to allow us an analytical treatment. We hud a strong dependence of the excitation spectrum 
upon the shape of the potential on the true vacuum side. We then discuss features of the excitation 
spectrum common to general tunneling potentials not restricted to our simple model. 
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There is a growing interest in the quantum state of a scalar held inside a vacuum bubble that appears after false 
^T^acuum decay. This is because the spacetime region inside the bubble may be considered as a homogeneous and 
j^isotropic open universe and, when combined with a class of inhationary models, there is a possibility that our universe 
entirely contained in a single bubble and the present large scale structure of the universe with a low density 
parameter, SIq ~ 0.1, may be explained [1-6]. In the standard scenario, inhation solves both the horizon and hatness 
problems simultaneously. But in the new one-bubble scenario, these problems are solved by two different inhationary 
stages of the universe. In the hrst inhationary stage, the universe is in a false vacuum. If this stage lasts long enough, 
the universe may be approximated by a de Sitter spacetime with high accuracy when the false vacuum decay occurs 
through quantum tunneling. This implies the decay is mediated by the 0(4)-symmetric Euclidean bounce, which 
consequently ensures that the spacetime after the decay has 0(3, l)-symmetry [7]. In other words, the region inside 
a nucleated bubble is a homegeneous and isotropic open universe. Then if the vacuum energy inside the bubble is 
nonvanishing, the second inflationary stage follows and a large amount of entropy will be produced at the end of the 
second inflation to solve the flatness problem.^ 


*JSPS junior fellow 

^Here by the flatness of the universe we mean Qo is not greatly different from unity. 



In order for this one-bubble inflation scenario to be successful, it is then a matter of great importance that if the 
quantum fluctuations of a scalar held inside a nucleated bubble can account for the observed large scale structures of 
the universe as well as for the detected anisotropies of the cosmic microwave background. As has been calculated in 
Refs. [4,5], the information of the quantum state at the beginning of the second inflation inside a nucleated bubble 
remains until today and is reflected in the CMB anisotropies on large angular scales, particularly in those corresponding 
to supercurvature scales at matter-radiation decoupling. It is thus of particular interest to clarify properties of the 
quantum state inside a nucleated bubble and how they are brought forth through false vacuum decay. 

In one of our previous papers [15], we have investigated this problem by assuming that the quantum state is in the 
Euclidean vacuum (Bunch-Davies vacuum). This means we have neglected the effect of quantum excitations through 
tunneling process but only taken into account the effect of the background spacetime curvature. Although one can 
construct a model in which this is a good approximation, such as the one advocated recently by Linde [6], it will not 
be so if there is only one scalar held and it is responsible both for the false vacuum decay and for the second inflation. 

On the other hand, efforts to understand the quantum excitations during tunneling process was Rrst made by 
Rubakov [8] and then by Vachaspati and Vilenkin [9]. Subsequently, based on the multidimensional wave function 
approach [10,11,14], we have investigated quantum excitations produced through false vacuum decay [12,13]. In 
particular, in [13], we have presented a systematic method to evaluate the degree of quantum excitations and calculated 
the excitation spectrum for simple models. However, there we have made several assumptions which may not be 
relevant for the one-bubble inflation scenario. One of them is the neglection of the background spacetime curvature. 
Another is that we did not consider the excitations of the tunneling scalar held itself but some other one that couples 
to the tunneling held through the mass term. Although these simplifications made it easy to understand gross features 
of the quantum state inside the bubble, we have to admit that it is far from complete. For example, if one considers 
the excitations of the tunneling held, there exists inevitably a region of negative mass-squared around the top of the 
potential barrier and that may affect the results seriously. 

In this paper, as a step toward full understanding of the matter, we extend our previous analysis [13] and investigate 
the quantum excitations of the tunneling held itself through false vacuum decay. However, for simplicity, we neglect 
the background curvature and assumes the Minkowski background. The paper is organized as follows. In section 2, to 
make the paper self-contained, we briefly review our formalism for determining the quantum state after false vacuum 
decay. Our formalism uses the fact that the excitations can be concisely described in the language of particle creation 
due to a varying mass, as that in curved spacetime. In section 3, we present a model of the tunneling potential which is 
piece-wise quadratic, hence is simple enough to allow us an analytical treatment but is expected to retain the essential 
feature specific to the excitations of a tunneling held. Then we express the resulting quantum state in the language of 
particle creation and give a complete (but complicated) formula of the particle spectrum for this model. In section 4, 
based on the formulae derived in section 3, we evaluate the particle spectrum for various model parameters in detail. 
For limiting values of the parameters, we give analytical expressions for the particle spectrum. For other values of the 
parameters we show the results obtained numerically. In section 5, we discuss possible role of discrete modes which are 
associated with the vibration of the bubble wall. Although these modes cannot be interpreted as particle modes, we 
argue that they will also contribute to the quantum state inside the bubble but those with spherical harmonic indices 
f = 0 and 1, which represent translational degrees of freedom of the bubble location and gives rise to divergence in 
the two-point function, will be absorbed into the metric perturbation when gravitational degrees of freedom are taken 
into account. Finally, section 6 is devoted to summary and future issues. 

Throughout this paper, the metric signature for Lorentzian spacetime is (-h -f-f) and the units h = I and c = 1 

are used. 


II. REVIEW OF FORMALISM 


To begin with, we give a brief sketch of our method for solving the evolution of the quantum state through false 
vacuum decay. We use the multidimensional WKB wave function approach originally developed by Banks, Bender 
and Wu [10] and de Vega, Gervais and Sakita [11], and recently elaborated by us [14]. We consider a scalar held (f) in 
the Minkowski spacetime whose action is given by 


S' = 





( 2 . 1 ) 


where V{(j)) is a tilted double-well type potential as shown in Fig. 2. We consider the situation in which the held is 
initially at the false vacuum {4> = 4>p) and decays toward the true vacuum {4> = 4>t) by quantum tunneling. 

In order to trace the evolution of the quantum state through false vacuum decay, we construct a quasi-ground state 
wave functional which is an energy eigenstate of the time-independent Schrddinger equation and is exponentially close 


2 



FIG. 1. A schematic picture of the bounce solution. 
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In the next WKB order, quantum fluctuations around the bounce solution come into play. Setting (f) = (f)B + the 
WKB wave functional to this order is given by [11,14] 


'?[</'(•)] =W) exp 


dT' / 


2 I (9T' / 2 I (9x 




exp 


J J d^xd^yip(x)n(x,y;T)ip(y) 


(2.5) 


where 


fi(x, y; T) := J , T) . 


( 2 . 6 ) 


The function g\^(x,T), which we call the mode function, satisRes the field equation for Lp on the background <f>B, 


i = l 


fifk(x,T) = 0 


(2.7) 


and the inverse g-^ is defined by 


/ d^x^k(x,T)^j^/(x,T) = - k') 


( 2 . 8 ) 


To solve Eq. (2.7), we need to set an appropriate boundary condition. At false vacuum 4>b = ff'F (T = t ^ ~co), the 
second line of the wave functional (2.5) must coincide with that of the would-be ground state at false vacuum. This 
requirement determines the boundary condition at r ^ —oo as 

(,k(x,r)^e-»‘^’^+“-'V (2.9) 

where Wk := + V {<f>F) ■ In the Lorentzian region A4, the function gk{x, it) is given by the analytic continuation 

of (/k(x, r) through r = 0. After all, the problem of constructing the WKB wave functional reduces to the problem 
of solving Eq.(2.4) and Eq.(2.7). 

In A4, the state of Lp described by the second line of Eq.(2.5) has a simple interpretation. In the second quantization 
picture (regarding t as real physical time), p is represented as a field operator, 

<,c>(x,t) = y d^k yk(x,t)&(k)-b Mk(x,t)&^(k)^ , (2-10) 

where &k and are the annihilation and creation operators, respectively, relative to a state |0) annihilated by &k, 
and a bar denotes complex conjugation.. The mode functions Mk(x,t) satisfy the field equation and are normalized 
with respect to the Klein-Gordon inner product, but are not necessarily positive frequency functions. Hence the 
‘vacuum’ |0) annihilated by &(k) is not generally an eigenstate of the Hamiltonian unlike the Minkowski vacuum |0). 
The former is a squeezed state over the latter, given by a Bogoliubov transformation, and contains non-vanishing 
spectrum of excited particles. It is known that the wave functional for the state |0) is given by 

(^(•)|o) cx exp 


j j 


( 2 . 11 ) 


where 


Q(x,y;f) := J d^k^^(y, f) . (2.12) 

This expression is nothing but the second line of Eq.(2.5) if Mk(x,t) is identified with g\^(x,it). Since Sl(x,y;t) is 
invariant under linear transformations of g\^(x,it), we then conclude Mk(x,t) = V)k' Ckk'dk'(x, *f) for some Ckk' with 
det Ckk' 7^ 0. Thus the quantum state of p after tunneling is a squeezed state over the Minkowski vacuum determined 
by the mode function Mk(x,t). 
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In order to find the mode function Mk(x,f), it is convenient to use the coordinates which respect the symmetry of 
the background <f>B, i-e., 0(4) in £ and 0(3, 1) in A4 [13]. In the Euclidean region £, we use , Xe, 0, v) where (6, cp) 
are the usual two-dimensional spherical coordinates, and (Ce,Xe) are related to r = |x| and r as 

r = ^EsmxE, r=-^EC0SXE, 

0<XE<^, 0<e£<oo . (2.13) 

Then the Euclidean metric takes the form, 

ds^E = ‘iCs + &{dXE + XEd£l^). (2-14) 

In the Lorentzian region A4, we have two distinctive spacetime regions characterized by the action of 0(3, 1). They 
are separated by the light cone that expands from r = t = 0. The spacelike region is called the (spherical) Rindler 
space, and the (future) timelike region the Milne universe. Since we are interested in the quantum state inside the 
bubble, we focus on the region corresponding to the Milne universe. Then the coordinates in the Milne universe are 
obtained by the replacement (Ce,Xe) (~*^)*x)) which yields 


r = ^ sinh X, i C cosh x, 

0<^<oo, 0<x<oo . 

The metric of the Milne universe is given by 

dslf = + i^{dx^ + sinh^ xd£i^)- 


(2.15) 


(2.16) 


It is to be reminded that 4>b is constant on the = const, (or ^ = con,st.) hypersurface. With these coordinates the 
Laplacian (or d’alambertian) which appears in the held equation is rewritten as 


dr^ 


E 

i = l 
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52 
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9**2 d^l ^E dis i 


2 ^E ) 


T 2 .— 

LiE — 
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sin Xe ^Xe 


sm XE 


dXE 


sm XE 


T ^ 


(2.17) 


in the Euclidean region, and as 

32 3 


9f2 


E 


92 


3 9 
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9**2 


L2 := 


de e 
1 9 /^. ,2 9 

smh X ox \ 9x 


• 1 2 ’ 

Sinn X 


(2.18) 


in the Lorentzian region, where is the Laplacian on the unit two sphere. 
As the 0(4)-symmetric bounce solution is considered, Eq.(2.4) reduces to 


r d2 

deE 


_3_ d 

^E d^E J 


f>Bi^E) = 14'(9_b(^_e)) 


in £. The boundary condition is 


t>Bi^E) <f>E for ^E 00 , 


(2.19) 


( 2 . 20 ) 


In A4, the equation for 4>b is obtained by replacing by —in Eq.(2.19), which describes the evolution of 4>b inside 
a nucleated bubble. 

Given the background solution <f>B, we need to solve for g\^. Eor this purpose, we expand g\^(x,it) in terms of 
harmonic functions on the three dimensional unit hyperboloid. 
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YplmiXj ; fplix) — 

which is an eigenfunction of the Laplacian operator — 

_T 2 -y _ 

^ plm — 


T(ip + £+l) 


T(ip) 


1 p-^-l/2 


with eigenvalue 1 + 


(1 + p'^)Yp(„ 


(coshx), 


Then Pk(x, it) is expanded as 


fi'k(x, it) 


roo _^ _ _ jp _ 

I dp'^X{k;pim)upi,j,{^,xMy, UplmiLxM) = ^p Ypimix, 

° im ^ 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


where we have associated the complex conjugate of Up£m with pt in accordance with the general discussion. We 
mention that p corresponds to a comoving wavenumber in the Milne universe and p = 1 to the spatial curvature scale 
on the ^ = const, hypersurface. We then consider the analytic continuation of Upjm to the Euclidean region by ^ ^ i^E 
(0 < ^e) and x ^ ~iXE (0 < Xe < "^/S)- In £, Eq.(2.7) for pt reduces to the equation for Gp(^E) '■= Fp(i^E), 


r 52 
Ml 


Id p2 

^E d^E il 


y" {dBiiE)) 


GpiM) = 0 


(2.24) 


It has been shown that the boundary condition (2.9) for g\^(x,T) at r ^ —00 means the regularity of pt on the 
r < 0 half of the complex T-plane. Since fpi{—iXE) is regular for 0 < xe < 7r/2, this condition is translated to the 
boundary condition for Gp(^E) that [13] 


GpM) Kipim^E) for ^e 00 . 


(2.25) 


In X4, Fp(^) is given by solving Eq.(2.24) with ^e replaced by —i^ and with the asymptotic boundary condition 
Fp(^) = Gp(—i^) at ^ ^ 0. As ^ ^ 00 , the bounce solution dsiO undergoes dampted oscillations around the true 
vacuum d = dr- Therefore at ^ ^ 00 , Ep(^) will generally have the form, 

PpiO = —r^MfpiO + —^'MfpiO ’ (2.26) 


where 


fp(0 = ^GPGHG\MO; M^ = V"{dT). (2.27) 

Note that fpiP) is the positive frequency function for the Minkowski vacuum [13]. That is, one can expand the held 
operator as 


fiLxM) 


dpY, 


fpiO 

c 


Ypi m {X,^)ap£ 

m H“ 


fpiO 

c 


YpimiXj 


(2.28) 


and the Minkowski vacuum is annihilated by apim- On the other hand, as mentioned previously, the quantum state 
after tunneling is a “vacuum” state |0) specihed by regarding Fp(d) as the (unnormalized) “positive frequency” 
functions. The orthonormalized positive frequency mode functions are then given by Upim with Fp replaced by 


UpiO 


4 (|ci,| 2 -|c 2 ,| 2 )^"(^^- 


(2.29) 


Thus the spectrum of created particles Up observed by the Minkowski vacuum observer in the asymptotically future 
region is given by 


= npb{p-p'l 


Clp/C2p 1 


(2.30) 
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FIG. 2. The potential of the tunneling held. 

Here we need to comment on the scaling property of this system. Under the rescaling given by 

^ gcfii (i = l,2) 

X g'x, 

m^m/rj, n^n/g', M^M/g', 


(3.4) 
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the action scales as 


2 


(3.5) 


S' 


g' 


s. 


Thus the system is transformed to the same system but with a different Planck constant. Since the decay rate of the 
false vacuum is determined by the value of the action of the bounce solution, the tunneling rate changes under this 
rescaling. However as the held equation is invariant, the particle creation rate does not change. Thus every rescaled 
model is equivalent with each other as long as our interest is restricted to the particle creation. Therefore specifying 
the only three non-dimensional parameters, say m/M, ji/M and 4>2/4>i is sufhcient for our purpose. 


A. bounce solution 


First let us consider the bounce solution. In our potential model, we may regard the bubble wall to be the spacetime 

region surrounded by the spheres at and ^2 where (pi = (pBiCi) (* = 1,2). So let us call and ,^2 the inner and 

outer radii of the wall, respectively. With this dehnition, the wall is nothing but the negative mass-squared region. 
Here we tacitly assumed that the wall is entirely contained in the Euclidean region 5; 0 < < ,^2 < co. But there 

exists a case in which the wall extends to the Lorentzian region A4. Since the scalar held (pB can get over the potential 
barrier only in £, the outer edge of the wall ,^2 must be in 5, but the inner edge of the wall can be either in £ or 

in A4. We call the former the EE case, while the latter the EL case. We Rrst consider the EE case. The EL case will 

be discussed at the end of this subsection. 

For the potential (3.1), Eq.(2.19) reduces to Bessel’s differential equation. Therefore the bounce solution 
which satisfies the boundary condition (2.20) takes the form. 


<?^1 


, Bi{}i^e)~^ hi^^E) + B2 {}i^e)~^N iiji^E) + — 6 < ■Sb < 6 , 

A* 

C{miE)~^Ki{miE) -6 < ■Sb < 00 , 

V r^Z 


(3.6) 


where and are the Bessel functions of the first and second kinds, respectively, and and are the modified 
Bessel functions of the first and second kinds, respectively. There are six unknown (non-dimensional) variables in the 
above; A, Bi, B 2 , C, M^i and M^ 2 - 

As the potential is constructed to be smooth to its first derivative, (pB, d(pB/d^E and d'^<pB/d^'^ must be continuous 
everywhere. Requiring this continuity at the junction points ^e = Ci, together with the equalities (pi = <?^b(6 ') 
(i = 1,2), leads to the following algebraic equations: 


1 = A(M6)-'7i(M6) = (//6)“' [BiJM) + S2Ai(//6)] + 4 > 

A* 

AhiM^i) = -Bi J2(//6) - B2N2{gi^i), 

Ah{M^,) = 4 [BiMb^i) + S2AV3(a^6)] , (3.7) 


and 


-^ = C{m^2) ^Ki{m^2) -= {l^i2) M-8i7i(/^6) +-B2AVi(//6)] + 4 ’ 

Cj) ^ TT~t j-I) 

CA2(m6) = Si J2(//6) + S2A2(/«6), 

CA3(m6) = £ [Si J3(//6) + S2A3(/«6)] , (3.8) 

where the second equality in the first line of each set of the equations is the continuity of (pB, the second line is that 
of d(pB/d^E, and the third line is that of d/d^Ei^E^d(pB/d^E)- Of course, these equations are not independent with 
each other; the third lines of each set can be derived from the rest of equations. Hence there are six independent 
equations which are necessary and sufhcient to determine A, Si, S 2 , C, M^i and M^ 2 - 

The analytic continuation of this bounce solution to A4 is given by the replacement of ^e by —i^ as before: 

<PB(-iO=A(MO-^Ji(MO<Pi . (3.9) 
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As the oscillation around the true vacuum attenuates, 4>b does not reach the junction point 4>i any more in A4. Hence 
(^B is conhned in the region V"((^) = and there is no additional particle creation after tunneling in this case. 

Here it is to be noted that the equation for the mode function Gp(^E) (2.24) is expressed not in terms of the 
original potential parameters m/M, ji/M and 4>2/4>i but rather of parameters which specify the wall conhguration. 
In the present case, the wall conhguration parameters are m/M, fi/M, and M^ 2 - Thus instead of the original 

potential parameters we may regard these four parameters as those specify the model, three of which can be given 
freely. A convenient choice for the discussion given below is to give m/M, M^i, M^ 2 - With this choice we hud ji/M 
is determined from the four equations consisted of the second and third lines of Eqs. (3.7) and (3.8), despite the fact 
that there are hve unknowns; A, Bi, B 2 , C and fi/M. This is because these four equations are homogeneous with 
respect to A, Bi, B 2 and C. In fact, from the condition that there exists a non-trivial solution for A, Bi, B 2 and C, 
we can derive the equation to determine n/M, 


det 


( /3(M6)J2(//6) + ^/2(M6)J3(//6) /3(M6)#2(//6) + ^/2(M6)#3(//6) 

Ks{m^2)J2{lJ-^2) -if2(m^2)=/3(/u6) A'3(m^2)772(/u6)- K2{m^2)Ns{iJ.^2) ' ~ 

\ m m 


(3.10) 


The fact that n/M is determined in this manner can be explained by taking the derivative of the equation for the 
bounce (2.19), 


I -lT/'Vc i 


i’ = o, 


where t/> := d(f)B / dS,E- The function t/> satisRes the boundary condition, 

t/i ^ 0 for ^ 0 and 00 , 


(3.11) 


(3.12) 


and it has no node. This is just the condition that the ground state wave function for a one-dimensional quantum 
mechanics with the potential l'b/{A(//^)-\-V"{(,E) should have zero energy eigenvalue. Given m/M, and M^ 2 , this 

condition is exactly what is expressed by Eq.(3.10). Moreover, from the above analogy with the quantum mechanics 
it is clear that /j, cannot exceed the order of 1/(^2 ~ ■fi)- 

Now let us discuss the EL case in which the inner edge of the wall is in the Lorentzian region A4. In this case the 
bounce solution (/b is given by 


(/Bi^E) _ f + ^; 0<£.e <£.2, 

[CimCEy^KiimCE) -; £2 <£e<oo, 


(3.13) 


with the junction conditions. 


-^ = C{m£2) ^ifi(m6) -= (/^6) ^BJi{h£2) + , 

(pi m^ 

CK2{m£2) = 5J2(/u6), 

CAG(m6) = -BUji£2)- (3.14) 

m 

By the same reason as in the EE case, we have three independent equations and there are the same number of 
unknowns; B,C and M£ 2 - 

The analytic continuation of this bounce solution to A4 is 

_ i ; 0<^<6, 

-^1 \ai(m0-^MM0 + A2(m0-^Ni(m0; £i<£<co, 


with the junction conditions, 

1 = B(fi£i)-^Ii(fi£i) + 4 = Ai(M£i)-^Ji(M£i) + A2(M£i)-^ Ni(M£i), 

A* 

Bl2i^£i) = AiJ2{M£i) + A2N2{M£i), 

M 

Bhi^i£i) = — iAiJ3{M£i) + A2N3 {M£i)) . (3.16) 

A* 
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Here again only three of the above equations are independent and the unknown parameters are Ai, A 2 and M^i, 
assuming B is known by solving Eqs.(3.14). 

As in the EE case, we may choose the wall conhguration parameters, m/M, M^i and M^ 2 , as independent model 
parameters. As before, ji/M is determined from the last two equations of (3.14). Specihcally the equation to be 
solved is 


A3(to^2)=/2(/u 6) = ^A'2(m^2)=/3(/u6)- (3-17) 

Thus in the EL case fi/m is independent of M^i and is given as a function of m^ 2 - More precisely, a close analysis of 
Eq. (3.17) reveals that 12^2 varies from ji i ~ 3.8132 to j 2 q ~ 5.1356 as m ^2 varies from zero to inhnity, where jm,n 
is the n-th zero point of the Bessel function Jm{z)- 


B. mode functions and particle spectrum 

As we have found the background solution <f>B, we now turn to the equation for the mode function (2.24). As in 
the case of the bounce solution, Eq. (2.24) reduces again to Bessel’s differential equation. Eirst consider the EE case. 
In this case, Gpl^s) in £ is solved to give 


( aipKip{M^E) + a.2pIip{M^E)', 0 < , 

Gpi^E) = < bipJip{n^E) + b2pNipln^E); < ^2 , (3.18) 

[ Kip(mCE); &2 <Ce < 00 , 


where the function Kip is chosen in the interval ^2 < ■Cb < co in accordance with the boundary condition (2.25). This 
time, Gp(^e) and its Rrst derivative are necessarily continuous at the junction points. These conditions determine 
the coefhcients aip, a 2 p, and & 2 p- Using the matrix notation, aip and a 2 p are expressed as 


f aip 

V «2p 


= XpYp 


Kip{m^2) 


where 


^ f /^ 77 /p(// 6 ) -NipifiS)] 

"" 2 J^p(A2) )■ 


(3.19) 


(3.20) 

(3.21) 


The mode function Fp(^) in the Lorentzian region A4 is given by Eq. (2.26), where the coefhcients cip and C 2 p are 
related to aip and a 2 p as follows. Eollowing the prescription ^ analytic continuation of the Bessel functions 
gives 


IGp{-iM0 = ye-"W2^g)(M0 , 
I,pi-iM0 = (AggMO + AggMO) 


This implies the following relation, 

( GP \ f flip A 

ycijiy 2v0 1 J \a2p J ' 

Since 

Gpi^E) = (aip - sinh( 7 rp)a 2 p^ Kip{M^E) + ^ [U>(Af^£;) + I-ip{M^E)] , 


(3.22) 


(3.23) 


(3.24) 


for 0 < where both Kip and Iip-\-I-ip are real, the fact that Gp(^E) is real in £ implies both aip — sinh(7rp)a2p 

and a 2 p must be real. Using this property, the expression (2.30) for the particle spectrum Up is rewritten as 
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(3.25) 


^ TT^jaipp 

where aip and a 2 p are given by Eq. (3.19). This expression, although exact, is too complicated that it is almost 
impossible to gain any physical insight from it. In the next section, we will examine several extreme situations in 
which we can obtain approximate analytical formulas which are more comprehensible. We will also evaluate Eq.(3.25) 
numerically to Rll the parameter regions not covered by the approximate formulas. 

Now we consider the EL case. In this case, the mode function Gpl^s) and its analytic continuation Fp(^) = Gp(—i^) 
are expressed as 


Gpi^s) 

FpiO 


Kipim^E)] ^2<^E<oo, 

bipJip(l^CE) + b'jpNipifi^E)', 0 < Ce < C'j , 

bipIipijJ-O + hpKipiM)'^ 0 < ^ < , 


(3.26) 

(3.27) 


Through the origin ^e =■5 = 0, the Bessel functions are analytically continued as 


M-iK) = , 

Nip{-iK) = -iG^''^Iip{K) - -e-^^''^Kip{K) • (3.28) 

TT 


With the aid of these relations, we get 


where 


^Ip 


\ 


— 7 W Y 

— p ^ p 


Kip{m^2) 

mKlp{mi2) 


Tip 

Wp 


TT^i / Me^PGH(%{M^i) -eGPGHG\MCi) 

-4i -Me-^PGH(%{M^i) e-^PGHF\M^i) 

f e^PG -ie^p/"^ \ 

V 0 -fe—p/2 J ’ 


Iip{p-£,l) Fip{jJ,S^i) 


(3.29) 


(3.30) 

(3.31) 


and Yp is given in Eq. (3.21). Substituting cip and C 2 p given by Eq. (3.29) into Eq. (2.30), we obtain the expression 
for the spectrum of the created particles in the EL case. 


IV. DETAILED ANALYSIS OF PARTICLE SPECTRUM 

As shown in [13], the spectrum of created particles has the general feature that rip is nearly constant for 0 < p ^ 1 
and decreases exponentially as e^^’^Pfor p ^ 1. Hence the particle spectrum is basically obtained if rio := np=o is 
known. Note that p = 1 corresponds to the curvature scale of the ^ = const, hypersurface in the Milne universe. It 
should be also noted p = 0 does not corresponds to the zero mode of the ^ = const, hypersurface but to the mode 
with characteristic scale of the curvature as well. Furthermore, if one considers implications of the present analysis 
to the one-bubble inflation scenario, what one wants to know most is the curvature perturbation spectrum on large 
scales comparable to the spatial curvature scale, which is described by rip at p < 1. Thus we may focus on the plateau 
of rip at p ^ 1. For definiteness, we take no as the representative value. 

For limiting cases in which the argument of the (modified) Bessel functions is very small or large compared to unity, 
we can derive a rather simple analytical expression for no by using the asymptotic behaviour of the (modified) Bessel 
functions which may be found in Ref. [16]. As discussed in the previous section, there are two ways of specifying pa¬ 
rameters of the model. The potential parameters are less directly related to rip than the wall configuration parameters. 
Therefore, in what follows, we first consider cases with extreme values of the wall configuration parameters. After 
disclosing the relationship between the wall configuration parameters and the particle spectrum, we then interpret the 
results in terms of the original potential parameters by analyzing the relations between these two sets of parameters. 
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FIG. 3. The typical shape of the particle spectrum np. 


A. no as a function of wall configuration parameters 


(i) Thin wall EE case: 

We first consider the case when fh^i and m ^2 ^ co, where fh represents either m, M or /r, and (^2 ~ ^ 1; 

which corresponds to the thin wall limit. We regard 1/m^i and l/m ^2 as small parameters, which we denote by e, 
and assume /u (^2 ~ ■fi) ~ 1- 

After tedious manipulations, the expansion of Eq. (3.10) with respect to e to 0(e) gives 

(Mm - ^?)T + + m) - ^ (— + ^ (m - /rT) + ^ ("— + —^ (M - ^T) = O (e^) , (4.1) 

8^1 \ ^ M J 8^2 \^i mj 


where 


T := tan/r(^2 - 6)- 

With the aid of this relation, the particle spectrum at p = 0 given by (3.25) is evaluated as 


no 


M M _ M 

p Mj \fiCi 


-2 


Equivalently, again using Eq.(4.1), no may be expressed as 


no 


M 


f ii{M + m) 

, arctan —r- 

M J \ _ Mm 


exp(—4M^i). 


(p^i)'^exp(-4M^i), 


(4.2) 


(4.3) 


(4.4) 


where arctan takes the value between 0 and tt. 

One sees that no is exponentially suppressed as the feature that has been found in the previous analysis of 

the thin wall limit [13]. However, one also Ends a large factor (p'fi)'^ in front, which were absent in the simple model 
discussed in [13]. 

(ii) Boundary between the EE and EL cases: 

The boundary between the EE and EL cases is given by the limit M^i ^ 0 and ^ 0. In this case the limit of 
either m ^2 ^ 0 or m ^2 ^ co can be treated analytically. 

As mentioned below Eq. (3.17), for m ^2 ^ 0, p is fixed as 


fi ^2 — ji,i ~ 3.8132. 


(4.5) 


Then we obtain 


no 


1 

4 


ki -\ -In 


M 


2 


12 


(4.6) 



where ki is a numerical factor given by 


ki= jii.iA^iO'i.i) (^In +'1^ +A^o(ii.i)| j JoO'i.i) 


-3.9245 In ( ^ 1 - 2.3929, 


and 7 ~ 0.5772 is the Euler constant. 
For m ^2 CO, /j, is hxed as 


Then we obtain 


/«6=i2.i~ 5.1356, 


no 


1 

4 


A ;2 + - In 

TT 


M 


2 


where ko is given by 


k'j 


No{j 2 ,l) 

Jo(i2,l) 


2.4602. 


(4.7) 


(4.8) 


(4.9) 


(4.10) 


Thus, in both of the limiting cases, no is of order unity unless the ratio M/ji or mjjjL becomes too large or too 
small. 

(iii) EL case with ^ oo: 

In the EL case with ^ oo and ^ oo, a simple expression for no can be obtained for m ^2 ^ oo or m ^2 ^ 0. 
First consider the limit m ^2 ^ 0. In this case /j, is given also by Eq. (4.5). After a straightforward calculation, we 
then Rnd 


no 


k'i + I + M'-^ 
Hr fiM 


exp(2//^i). 


(4.11) 


In the case of the limit m ^2 ^ oo, [i is given by Eq. (4.8). Then we obtain 


no 


kl + l I? + 
Hr fiM 


exp(2//^i). 


(4.12) 


Thus in both limits of m^ 2 , the particle spectrum has the exponential factor 
common factor for the EL case irrespective of the parameters. 


This suggests that it is a 


Summarizing the above results of analytical tractable cases, we expect that the gross dependence of no on the wall 
configuration parameters is 


no ~ exp(-4M^i) , 


for the EE case, and 


no ~ exp(2//^i) , 


(4.13) 


(4.14) 


for the EL case. 

In order to test our expectation mentioned above, we have numerically evaluated no for various values of the wall 
configuration parameters. Figures 4 and 5 show no/exp(—4M^i) for the EE case and no/exp(2//^i) for the EL case, 
respectively, as functions of the parameters and M^ 2 , for typical ratios of m/M; m/M = 0.1, 1 and 10. We see 
that the approximations as Eqs. (4.13) and (4.14) are better than one might have anticipated except for some special 
cases. 
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(c) 


FIG. 4. The contour plots of the particle creation rate rap=o versus and — Affi in the EE case. They are plotted 
for a few typical values of mjM, i.e., (a) mjM = 0.1, (b) mjM = 1 and (c) mjM = 10. 


14 



(c) 

FIG. 5. The contour plots of rap=o versus and M ^2 in the EL case. As before they are plotted for the same typical 

values of m/M. 

As for the EE case, we see a ditch in the hgure, but it is in some sense superhcial. As is also observed in Eig. 3, there 
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FIG. 6. 
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FIG. 7. The plot of the function f{x). 
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FIG. 8. 

If we release a particle from the same point 4 >i as in the original potential for the case of the potential V, it cannot 
reach the false vacuum (fip. In order to make it reach (fip, the particle must be released from a point cj) = (fij closer 
to the true vacuum, i.e., 4>i > (f)j. Then which is the value of xe at (fii, must become larger. Thus we conclude 
M^i > M^i for with Rxed m/M and ji/M. 

In contrast to the EL case, we were unable to find a simple function of the potential parameters which is directly 
related to the amplitude of particle creation. This is because the bounce solution in the EE case depends on the mass 
at the true vacuum M as well as m and fi. Therefore M^i is generally fully dependent on all the three potential 
parameters. However, it is still possible to relate the shape of the potential to the particle creation in a couple of 
limiting cases, which we will discuss below. 

As noted in the beginning of section 3, 4>2/4>i is bounded from above by the condition that V{4>f) > Eor 

4 >2/4 >i close to the maximum value, we have ^ 1, which is just the thin-wall limit. In this limit, the particle 

creation is exponentially suppressed as In other word, as the wall radius becomes larger, the quantum state 

inside a nucleated bubble becomes closer to the Minkowski vacuum state. In this case, the wall radius is known to be 
given by [7] 
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6 ~ 6 ~ 35i/AT/, 


(4.23) 


where 


and 




•4>t 


'F 


(4.24) 


AV" := V{<f>F) - V{<f>T). 


Thus the amplitude of particle creation no is approximately given by 


no ~ exp 


12M5i 

Ay 


(4.25) 


(4.26) 


On the other hand, 4 >2/4>i is bounded also from below in order for the potential barrier to exist. Now as 4 >2/4>i 
decreases, decreases and before the potential barrier vanishes becomes zero, which is just the boundary of the 
EE case and the EL case. Thus in this limit, the analysis of the ^ 0 limit of the EL case is also appropriate. 

Summarizing the above analyses for the EL and EE cases, we conclude that no is approximately determined in 
terms of the potential parameters as 


“ \exp(-12M5i/Ay); 

Before closing this section, one comment is in order. In this paper, we have considered a potential model with a 
constant mass around the true vacuum. However in order for our results to be valid, the mass M(^) on the true 
vacuum side is not necessarily strictly constant. The only restriction is that M(^) should vary sufhciently slowly. 
Namely if 


(4.27) 


jM(0 < M(0 , (4.28) 

the mode function evolves adiabatically and there will be no additional particle creation on the true vacuum side. 
Now in a simple version of the one-bubble inflation scenario, the tunneling held also plays the role of the inflaton held 
inside the bubble, and the mass of the inflaton held changes very slowly in the slow rolling phase. Thus our results 
are expected to give non-trivial implications to the one-bubble inflation scenario, at least the effect of tunneling to 
the spectrum of the inflaton held fluctuations is concerned. We plan to make a detailed investigation of this issue in 
a future publication. 


V. DISCRETE MODES 

So far, we have not carefully investigated the completeness and normalization of the mode functions for the descrip¬ 
tion of a quantum state of the held cp. In order to specify a quantum state completely, we need a set of all possible 
mode functions which have properly normalized Klein-Gordon inner products on a Cauchy surface. However, the 
spacetime inside the forward light cone of the center of a bubble, which is described by the Milne universe, does not 
contain any Cauchy surface of the whole Minkowski space. 

Thus we Rrst need to introduce new coordinates which cover the spacelike region outside the forward light cone and 
which respect the symmetry of the bounce solution. Such coordinates are known as the (spherical) Rindler coordinates 
and the spacetime covered by them is called the Rindler space. The metric of the Rindler space is given by 

ds"^ = d^"^ + (y—dxR + cosh^ ; (5-1) 

where the XR = const, hypersurface of the Rindler space is a Cauchy surface of the whole Minkowski space. These 
coordinates are related to the Milne coordinates as 

■fi? = , XR = X+Y- (^-2) 
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Now we have to evaluate the Klein-Gordon norms of the mode functions in the Rindler space. In 

order to do so, we need the analytic continuation of Upim to the Rindler space: 


'^plm — 


Gpi^n) 

Cr 


{XR)^im (^) 7 


(5.3) 


where 


VpiiXR) = fptiXR - *7r/2) 


r(ip + f + i)r(-ip + f +1) 


r(ip)r(—ip) 


1/2 


V-icoshxi? 




(5.4) 


We note that Vpi(xR) now plays the role of a positive frequency function for the state |0). In the same way as 
presented in Appendix A of Ref. [15], the Klein-Gordon inner product of Upim on the XR = const, hypersurface can 
be evaluated as 


poo 

{uplm, Rp'l'm') ■= i cosh'^ XR / cICrCr 

■Jo 


f du 


'plm 


V dxR 


’^p'l'm' '^plm 


dUp! i' fYi' 

dXR 




where 


= i cosh^ XR ( 


\dxR"^' dxR 




(5.5) 


(5.6) 


In the above, we have put p' = p when evaluating anticipating that the factor NpP should be zero if p yt p'. 


which we will show below. The evaluation of the factor is straightforward and we obtain 

w(i) 2p . , 

= — smh Trp. 

^ TT 


(5.7) 


For and p' > 0 (in fact we have p, p' > 0), the factor can be evalnated from the behavionr of Gp{£,R) near 

the origin alone. In the EE case, noting that Eq. (3.24) implies 


Gpiin) = - W-^!^7_,4Me«), 

smh wp smh wp 

and nsing the fact that i{cip — cosh7rpc2p) and C 2 p are real, we obtain 


(5.8) 


(5.9) 


Thus we have 


{'^pln 


1 Un'P'm.' ) — 


4 (|cip|2 - |c 2 p|^ 


-6{p - p)6wK 


(5.10) 


This agrees with the norm calculated on the ^ = const, hypersurface in the Milne universe, Eq. (2.29). We can show 
the same result holds also in the EL case. The reason why the Klein-Gordon inner product can be evaluated in the 
Milne universe, despite the fact that it does not contain any Gauchy surface, is due to the fast fall-off of the function 
Vpi{xR) with positive p^ at XR co- That is, the mode function Upim vanishes fast enough at future null inhnity of 
the Minkowski space. 

In contrast, any mode function with negative p^ has a divergent Klein-Gordon norm in the Milne universe. When 
the Klein-Gordon norm of a mode diverges, the normalized mode function vanishes. Hence it would not contribute 
to the quantum fluctuations of the held if it also vanishes fast enough at future null inhnity. However, for p^ < 0, the 
fall-off of a mode function at future null inhnity is not fast enough. Thus, we cannot claim that the Klein-Gordon 
norm evaluated in the Milne universe is equivalent to that evaluated on a Gauchy surface. Therefore it sometimes 
occurs that some of the modes with negative p^ do contribute to the quantum huctuations of the held. 

Specihcally, as seen from Eq. (5.8), Gp(^ij) with p^ < 0 has the behavior ~ <^pPr + near the origin = 0, 

where ap and fJp are determined by solving Eq. (2.24) with the boundary condition (2.25). This boundary condition 
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is a necessary condition for the normalizability also for the < 0 modes but not a sufhcient one. For a mode to be 
normalizable, there is an additional condition that jSp should vanish. For almost all the modes with negative jSp 
do not vanish. Hence they do not contribute to the quantum fluctuations. However, for certain discrete sets of p, l3p 
may happen to become zero. Then those modes become normalizable. Apparently, they are bound state modes and 
cannot be described in the particle picture. 

Nevertheless, we can evaluate the contribution of these modes to the two-point function (Wightman function) which 
characterizes the quantum fluctuations of the held. Let us denote these discrete values of imaginary p by Since 
the Wightman function is given by the summation of products of the mode functions, it can be divided into two parts 
as 


G~^(x, x') = GJ(*, x') + GJ(*, x'). (5-11) 

The two parts GJ and GJ are the contributions from the discrete modes and the continuous modes, respectively. 
They are given by 


G^{x,x') = ^ 


2 pn sinh7rp„W 


( 2 ) 


= —y 

97r2 


sinhKIC 1 Gp„(e,j)Gp„(e),) 


27r2 47r sin 7r|p„ | 


where 


#(2) — 


d^R 

Cr 


(GpAuyf 


and 


poo _ 

Gy{x,x')= / dp'^ 

Im 


4(|cipP-|c2,P) 


'^plm i^x'^Upi 

m ) 


1 


psinpC Gp{^R)Gp{^'pp) 


Jo 87r(|cip|2 - |c 2 pp) sinhC CrCr 

where ( is the Lorentz invariant distance between x and x', 


* _ f cosh X cosh x' — sinh x sinh x' cos 0 (in Milne universe), 

( — sinh xr sinh x'r + cosh xr cosh x'r cos 0 (in Rindler space), 

cos 0 = cos 6 cos 6 ' + sin 6 sin 6 ' coslp — 4>'). 


(5.12) 


(5.13) 


(5.14) 


(5.15) 


The second equalities of both Eqs. (5.12) and (5.14) are obtained after summation over I and m, and they manifestly 
show the 0(3, 1) invariance. 

In the present model, there exists at lease one series of such discrete modes. They are related to the perturbations 
of the wall location. The radial part of these mode functions is given by the derivative of the bounce solution as 

G^alliU) = Jr^^^ ■ (5.16) 

d^R 

It is easy to show that this mode function satisRes Eq. (2.24) with = —4. If we put p„ = 2i in Eq. (5.12) it 
diverges. However the divergence comes from the monopole {I = 0) and dipole {I = 1) parts of these modes, which 
just represent translations of the origin of the coordinates (so-called zero modes). In fact it is easy to check that these 
modes are represented by linear combinations of d<j)B/dx^ where are the usual Minkowski coordinates. Thus they 
should be removed. Then the remaining part of the Wightman function becomes finite but at the expense of losing 
the Lorentz invariance. Using the explicit form of the associated Legendre functions with special values of the indices. 


^i.-Y/2(cOshz) 


2 sinh Rz 
TT sinh z R 




2 R cosh Rz — sinh rz coth z 
TT sinh z r^r"^ — 1) 


(5.17) 
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the contribution to the Wightman function from these discrete modes is expressed as 


where 






cosh2C 11 , ^ , 8 . , . , , ' 

C -^ cosh C + - smh xr smh x r 


^wlll = / d^R^R4>%{^R). 


(5.18) 


(5.19) 


The last term in the square brackets is not Lorentz invariant. However, since the removed monopole and dipole modes 
simply represent the global translation of the bubble, this apparent violation of the Lorentz invariance should not 
affect observable quantities. Focusing on the wall fluctuations in the Rindler space, this point has been discussed in 
detail by Garriga and Vilenkin [17]. 

In the Milne universe, one can interpret the fluctuation of the scalar held Lp as the perturbation of the intrinsic 
curvature of the (f) = const, hypersurface. The scalar-type curvature perturbation is described by a single potential 
function TZ as [18] 


1 


- (v'v,- - + p; (3 - v'^v,) 


-(1 + p-^)TZY; + -(4 + 

o 


TZ(x) 


(5.20) 


where Vj- denotes the covariant derivative with respect to the metric on the unit 3-hyperboloid, and the second line 
is the harmonic expansion of the Rrst line with Y being the abbreviation for Ypim(Xj^) ^nd being the traceless 
tensor harmonics defined by 


y® =- 

j 1-Pp2 |_ 


1 


V®V,' - 


y =: 


1 -by 


V]Y. 


2 3 


We note that W® satisfies the equations. 


^^3 = -\^2^3y , = -(7 + t?)y; 


(5.21) 


(5.22) 


By an infinitesimal coordinate transformation ^ ^ ^ ^ T(x), we have 

ip{x) = ip{x) — ^bT{x), TZ{x) = TZ{x) — —T{x), (5.23) 

C 

on the new hypersurface [18]. Hence setting T = Lpl<f>B, we obtain 

n{x) = —^Lp{x), (5.24) 

K 9 b 

as the curvature perturbation of the (f) = const, hypersurface. 

Now for = —4, one finds the trace of the curvature perturbation vanishes. Furthermore, from the first equation of 
(5.22), we see that YJ becomes transverse. Thus the p^ = —4 scalar-type curvature perturbation happens to become 
transverse-traceless. This suggests that 7 = 0 and 1 modes do not contribute to the curvature perturbation and 
the term violating the Lorentz invariance in Eq.(5.18) will disappear when the two-point function of the curvature 
perturbation is considered. In fact, this can be explicitly demonstrated by operating with T>ij on it. Consequently, 
we obtain 


(^(3)i?®(*)^(3)i?f(*'))y = -4 = 


1 




-viv, 

! 3 ‘ 


,cosh2C 11 , . 

C—1-^ - ^coshC 

smh (, 3 


(5.25) 


which is manifestly Lorentz covariant. 

From the above discussion, it is anticipated that the quantum fluctuations of the scalar held inside the bubble will 
be better understood when we include degrees of freedom of the metric perturbation. This point will be disclosed in 
the future work. 
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Now let us examine if there exist discrete modes other than p = 2i. Since the mode (5.16) has no node, it is the 
eigenfunction for Eq. (2.24) with the lowest eigenvalue, i.e., the smallest . Hence other possible discrete modes must 
have the eigenvalues in the range —4 < p^ < 0 and have at least one or more nodes. Now in the thin-wall limit, it can 
be explicitly shown that such modes do not exist. Also in this limit, one Ends the mode function Gp with p = 0 has 
one node and diverges to minus infinity as ^ 0. Then as we vary the model parameters continuously, the value of 
the mode function Gp=o at = 0 should cross zero if there should appear a bound state mode with one node. This 
implies the divergence of no for a certain set of the model parameters. However, as seen from Figs. 4 and 5, we have 
found no divergence of no- Hence in our model, we conclude that there exist no additional discrete modes other than 
those related to the wall fluctuations. 

We will present a detailed analysis of these discrete modes in future publication. 


VI. CONCLUSION 

In this paper, we have investigated the self-excitation of a scalar held in process of its decay from a false vacuum. For 
this purpose we have considered a model potential which is piece-wise quadratic, hence allows us analytical treatments. 
We have interpreted the resulting quantum state inside a nucleated bubble in the particle creation picture. Then we 
have found the following features of particle creation. 

When the spacetime region of the bubble wall, which is defined as the region in which the mass-squared of the scalar 
held is negative, is confined in the Euclidean region £, the number of created particles per each mode is exponentially 
suppressed in the thin-wall limit and at most of order unity unless the mass scale at the true vacuum is exponentially 
small compared with that at the top of the potential barrier. On the other hand when the wall region extends to the 
Lorentzian region A4, i.e., the mass-squared at the center of the bubble is still negative, the particle creation can be 
significantly enhanced. In this case, we have derived an approximate formula (4.22) for the particle spectrum as a 
function of the model parameters which determine the shape of the potential, where is defined by Eq.(4.15). 

In addition, we have also considered the effect of a set of discrete modes which describe the oscillation of the bubble 
wall. From a careful analysis of the Klein-Gordon norms, we have argued that these discrete modes do contribute to 
the quantum state inside the bubble, though they cannot be interpreted as usual particle modes. However, since the 
monopole and dipole components of these modes corresponds to the spacetime translation of the bubble center, we 
have argued that a consistent treatment of these modes requires the inclusion of gravitational degrees of freedom into 
the analysis simultaneously. 

In view of the above considerations, the next step to be taken is to take gravity into account in the background bounce 
solution. A framework in this direction has been already done in Ref. [19]. Hence it should be fairly straightforward 
to extend the present analysis to the one on the curved spacetime background. Then the second step is to study the 
effect of gravitational degrees of freedom on the excitation of a tunneling scalar held. In this respect, we expect that 
a formalism developed in Ref. [20] for dealing with the fluctuations around the bounce solution with gravity should 
be useful. After this second step, the role of the discrete modes will be clearly and unambiguously understood. 

Once these issues are settled, we will be able to talk about the quantum state of the scalar held inside the bubble 
with confidence. In connection with the one-bubble open inflation scenario, we will be able to discuss quantitatively 
the influence of the quantum fluctuations induced by tunneling on the primordial density perturbations and on the 
CMB anisotropies on large angular scales. 
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